Phase diagram of the Gaussian-core model 
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We trace with unprecedented numerical accuracy the phase diagram of the Gaussian-core model, a 
classical system of point particles interacting via a Gaussian-shaped, purely repulsive potential. This 
model, which provides a reliable qualitative description of the thermal behavior of interpenetrable 
globular polymers, is known to exhibit a polymorphic FCC-BCC transition at low densities and 
reentrant melting at high densities. Extensive Monte Carlo simulations, carried out in conjunction 
with accurate calculations of the solid free energies, lead to a thermodynamic scenario that is 
partially modified with respect to previous knowledge. In particular, we find that: i) the fluid-BCC- 
FCC triple-point temperature is about one third of the maximum freezing temperature; ii) upon 
isothermal compression, the model exhibits a fluid-BCC-FCC-BCC-fluid sequence of phases in a 
narrow range of temperatures just above the triple point. We discuss these results in relation to the 
behavior of star-polymer solutions and of other softly repulsive systems. 
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It is common knowledge that crystallization is induced 
by the strong Pauli repulsion between inner-shell elec- 
trons, causing the effective interatomic potential to blow 
up at short distances. However, the existence of a ther- 
modynamically stable solid phase does not necessarily re- 
quire a singular repulsion for vanishing interatomic sep- 
arations. As a matter of fact, a finite square barrier, 
equal to a positive constant e for distances smaller than 
a given diameter a being zero otherwise, is an example 
of a bounded repulsion that supports a stable solid at all 
temperatures. 0, In this respect, a Gaussian-shaped 
potential 



v(r) = eexp(— r 2 /cr 2 ) , 



(1) 



is a more realistic finite-strength repulsion. The so-called 
Gaussian-core model (GCM) was originally introduced 
by Stillinger. [j| Such a potential, despite the fact that 
it is finite even at full overlap between the particles, is 
nonetheless perfectly admissible as an effective potential. 
Actually, it is used to represent the entropic repulsion 
between (the centers of mass of)self-avoiding polymer 
coils dispersed in a good solvent. |J,|5j,|6j Two distinctive 
features of the GCM, which are absent in the simpler 
penetrable-sphere system, are: i) the existence of a max- 
imum freezing temperature, T max ; ii) the occurrence, be- 
low T max , of reentrant melting into a dense fluid phase. 
Stillinger noted in his original paper that, in the limit 
of vanishing temperature and density, the GCM parti- 
cles practically behave as hard spheres with increasingly 
large diameter. In this limit, the fluid freezes into a face- 
centered-cubic (FCC) structure at a temperature Tf(p) 
that vanishes with the number density p (from now on, 
temperature and density will be given in reduced units, 



T* = k^T/e, where fee is Boltzmann's constant, and 
p* = per 3 ). Indeed, a straightforward calculation of the 
total energy of different cubic crystal structures shows 
that the FCC structure is favored, at zero temperature, 
only for reduced densities lower than 7r^ 3 / 2 ~ 0.1796. Be- 
yond this threshold, a body-centered-cubic (BCC) solid 
takes over. However, upon compression, any regular ar- 
rangement of particles is eventually destined to collapse 
for any T > (reentrant melting). |7j 

A comprehensive study of the phase diagram of the 
GCM was recently carried out by Lang and cowork- 
ers. 1^1 These Authors employed an approximate integral- 
equation theory to describe the disordered phase and a 
variationally-adjusted harmonic interaction for the crys- 
talline phases. The resulting phase diagram (see Fig. 9 of 
Ref. |^) accounts for the existence of a fluid phase and 
two solid phases. More specifically, the solid was found to 
be thermodynamically stable for temperatures lower than 
^max = 0.0102, a temperature at which the freezing line 
attains its maximum value for /9* nax = 0.2292. The phase 
diagram displays a fluid-BCC-FCC triple point that was 
estimated to fall at T*. = 0.00875. The FCC phase turns 
out to be stable below T t * in the low-density region (for 
p* less than « 0.17), whereas the BCC structure prevails 
for larger temperatures and densities. 

A parallel numerical study of the phase diagram of 
star-polymer solutions in a good solvent was carried out 
under a different assumption for the effective pair po- 
tential, modeled with an ultrasoft logarithmic repulsion 
within a diameter a and with a Yukawa potential outside 
the core. In this case, Monte Carlo (MC) simulations 
and free-energy calculations lead to a phase diagram that 
shows a rather complicate interplay between various cu- 
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TABLE I: Excess Helmholtz free energy per particle / cx , in 
units of k B T, calculated for some FCC (TV = 1372) and BCC 
(TV = 1458) solid states of the GCM. For T* = 0.003 and 
0.006, the tabulated values refer to systems with 864 and 1024 
particles, respectively. For each state and solid structure, we 
also display (within square brackets) the value of the reduced 
elastic constant c* = c<r 2 /e that intervenes in the Frenkel- 
Ladd calculation: for the given c, the mean square displace- 
ment of the Einstein crystal approximately matches the mean 
square deviation of a GCM particle from its reference crystal 
site. For a number of selected p* and T* , we verified that 
the quantity /3f«x.(N) + kxN/N, with f3 = (fc B T)"\ scales lin- 
early with TV -1 for large TV, in agreement with a conjecture 
formulated in pd|. 
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bic phases. It is desirable to have a similarly fully-Hedged 
analysis also for the GCM. 

To this purpose, we performed standard Metropolis 
MC simulations of the GCM, keeping the number of par- 
ticles TV, the volume V, and the temperature T con- 
stant. We used the particle-insertion method |i| and 
the Frenkel-Ladd technique [lflL ITlj to calculate the "ex- 
act" free energies of the dilute fluid and of the solid 
phases, respectively. By this means, we discovered that 
the phase diagram of the GCM is more elaborate than 
previously reported, showing elements of similarity 
with the phase diagram of star-polymer solutions. 

Our samples typically consisted of 1372 particles for 
the fluid and the FCC solid, and of 1458 particles for the 
BCC solid. Occasionally, we considered smaller as well as 
larger sizes, so as to check whether our conclusions were 
possibly undermined by a significant finite-size depen- 
dence. We payed much care to a safe estimate of statisti- 
cal errors. This is actually an important issue whenever 
different crystalline structures so closely compete, as in 
the present case, for thermodynamic stability. We com- 
puted free-energy differences between two equilibrium 
states of the system belonging to the same phase through 
standard thermodynamic integration. This method al- 
lows one to obtain the properties of the model for any 
state, provided that the absolute free energy has been 
autonomously computed in at least one reference state 
per phase. Table 1 gives the excess Helmholtz free en- 
ergy for some FCC and BCC states of the model. Though 
we did not systematically check the relative stability of 
other crystalline structures, we verified that, for tempera- 
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FIG. 1: Difference between the chemical potentials of 
pairs of GCM phases plotted as a function of the pressure 
along the isotherm T* = 0.002: pAfiauid, fcc (dashed line), 
/3A^fl uid ,BCC (dotted line), and (3A(i-pcc, bcc (continuous 
line). Upon increasing P (or p), the fluid transforms into 
a FCC solid; then, a FCC-BCC transition takes place until 
the BCC melts into a fluid phase again. The lines are spline 
interpolants of the data points. A zoom on the low-pressure 
region shows the non-monotonic behavior of /3A/iFcc, bcc, a 
feature that is ultimately responsible, at higher temperatures, 
for the reentrance of the BCC phase (see also Fig-EI - 



tures close to the triple point, the hexagonal-close-packed 
(HCP) solid is slightly less favored than the FCC solid, 
while the simple-cubic solid is not mechanically stable. 

The thermodynamically stable phase, for given tem- 
perature and pressure, is the one with lowest chemical 
potential fi(T,P). Figure^ shows the difference between 
the chemical potentials of competing phases plotted as 
a function of P at fixed temperature. The sequence of 
phase transitions undergone by the GCM at T* = 0.002, 
with increasing pressures, is fluid-FCC-BCC-fiuid: the 
disordered phase is actually seen to reenter the phase dia- 
gram at high density. As is also apparent from Fig.Q] the 
fi gap between the FCC and BCC phases is anything but 
monotonous when plotted as a function of the pressure. 
Note that, at low pressures, the BCC phase is about to 
become stable as the fluid freezes into a FCC structure. 
In fact, this eventually occurs at higher temperatures, in 
a way not documented before for the GCM. 

Upon increasing the temperature, the shallow valley 
in //fcc ~ Mbcc moves gradually upwards until a narrow 
range of pressure appears (for T* > 0.0030) where a sta- 
ble BCC phase slips in between the fluid and the FCC 
solid phases. This is possible because the corresponding 
increase of the freezing pressure with temperature is not 
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FIG. 2: Phase diagram of the GCM in the (p, T) plane, with 
a zoom on the triple-point region (inset). Transition densi- 
ties for each phase are shown for various temperatures: fluid 
(TV = 1372, black open circles), FCC (TV = 1372, blue solid 
squares), and BCC (TV = 1458, red solid triangles). Black 
crosses refer to simulations of larger systems (TV = 2048, fluid 
and FCC; TV = 2000, BCC) at T* = 0.0037. Continuous lines 
drawn through the data points are a guide for the eye. We 
also plotted the freezing and FCC-BCC coexistence loci cal- 
culated in |4J (black dashed lines), and the ordering threshold 
predicted by a one-phase entropy-based criterion (blue dotted 
line). 



large enough to suppress the reentrant BCC phase. For 
T* > 0.0038, the FCC phase ceases to be stable and a 
more regular behavior sets in, similar to that predicted 
by Lang and coworkers. [J] 

All in all, the phase diagram represented in Fig.|3 
emerges. If compared with Fig. 9 of |5], two differences 
stand out: a definitely lower triple-point temperature 
(T t * ~ 0.0031) and the as yet unpredicted reentrance of 
the BCC phase when the FCC solid is isothermically ex- 
panded for reduced temperatures in the 0.0031 — 0.0037 
range. In order to check whether this latter feature is 
a spurious effect due to the finite size of the system, 
we investigated the BCC-FCC phase coexistence also for 
larger samples, but we did not register any significant 
change in the location of the transition points (see the in- 
set of Fig.|2J|. In the triple-point region, the density jump 
is ~ 0.002 across the fluid-solid transition and ~ 0.00015 
across the solid-solid transition; the corresponding (abso- 
lute) values of the entropy jump per particle are ~ 0.7 fcs 
and ~ 0.1 respectively. The freezing line attains its 
maximum value (T^ ax ~ 0.00874) for /4 ax ~ 0.239. At 
the maximum, the fluid-solid transition is still first-order 
with an entropy gap between the two phases equal to 



FIG. 3: Temperature evolution of /3AiiFcc, bcc, the 
chemical-potential gap between the FCC (TV = 1372) and 
BCC (TV = 1458) phases, plotted as a function of the pres- 
sure. For T* — 0.003 and 0.006, the systems investigated were 
smaller (864 and 1024 particles, respectively). The dotted line 
refers to a calculation carried out at T* = 0.0037 using larger 
samples (FCC: TV = 2048; BCC: TV = 2000), and was plot- 
ted for a comparison with the smaller-sizes calculation at the 
same temperature. 



0.79 k B . 

To have a clue on the extent to which statistical errors 
may affect our conclusions, we turn the reader's atten- 
tion back to Fig.^ Let A/iA, b(-P) be the (i gap (at fixed 
temperature) between two generic phases A and B. In the 
triple-point region, we estimated a maximum statistical 
error on the minimum value of /3A/xfcc, bcc(-P) approx- 
imately equal to 10~ 2 . Of the same order is the maxi- 
mum error we estimated, near coexistence, on the values 
of /SA/Xfluid, solid However, the rate of change of this 
latter quantity is much larger, implying that its zero is 
more sharply defined. This means that fluid-solid coexis- 
tence is numerically better defined than solid-solid coex- 
istence. If we follow the evolution of /JAitpcc, bcc(-P) as 
a function of T (see Fig.[3J), we realize that, over the whole 
stability region of the reentrant BCC phase, this quantity 
takes values that are of the order of the estimated numer- 
ical errors. However, we can safely argue that the true 
errors are in fact much smaller since, otherwise, we would 
have hardly obtained the very smooth behavior repre- 
sented in Fig.|21as well as the clear phase portrait shown 
in the inset of Fig-El The already mentioned absence 
of any significant size-dependence of /3A/xfcc,bcc(-F > ) at 
T* = 0.0037 is a further guarantee of the reliability of 
the present phase diagram. 

The FCC phase of the GCM is energetically favored at 
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low densities (p* less than » 0.17), for temperatures up 
to T* = 0.008. This may actually explain why the FCC- 
BCC coexistence locus found in Ref. [j] is an almost verti- 
cal line in the p-T plane, which leads to a more extended 
region of FCC stability. In fact, Lang and coworkers used 
the Gibbs-Bogoliubov inequality to optimize a strictly 
harmonic model of both solid phases. This method may 
actually enhance the crystallinity and give, at the same 
time, an inadequate representation of the cntropic con- 
tribution to the solid free energies. Considering that the 
FCC-BCC transition occurs for rather small densities, 
the harmonic approximation seems a severe limitation of 
the theory. In fact, the variational technique basically 
propagates to higher temperatures the relative stability 
condition valid at T = 0. We also note that the FCC- 
to-BCC transition undergone, with increasing tempera- 
tures, by the GCM at low densities is to be ascribed to 
the higher entropy of the BCC phase, that is likely due 
to the presence of a larger number of soft shear modes. 
Even below the triple-point temperature, BCC-ordered 
grains tend to form in the liquid that is about to freeze, 
a phenomenon that substantially slows down crystalliza- 
tion. 

The low-density /low-temperature phase behavior of 
the GCM, with a triple point separating a region where 
the fluid freezes into a FCC structure from another region 
where these two phases are bridged by an intermediate 
BCC phase, is rather common among model systems with 
softly repulsive interactions, such as the inverse-power 
potential, v n (r) = A(a/r) n , flU IT3 and the Yukawa po- 
tential, v e (r) = Bcxp(-r/€)/r. E E3 The P hase 
diagram of the above two models is typically unfolded 
by one or two (possibly rescaled) thermodynamic quan- 
tities and by the relevant control parameter of the inter- 
action, i.e., the inverse-power-law exponent 1/n or the 
Yukawa length £. The crystalline pattern produced by 
such potentials is critically determined by their degree of 
softness, the BCC phase being promoted by a sufficiently 
soft interaction. A criterion to relate the phase behavior 
of v n (r) and vt(r) to that of the GCM is to require the 
logarithmic derivatives of such potentials to match that 
of the Gaussian potential (Eq.nj , at least for separations 
close to the mean interparticle distance, r — p~ x / z . The 
values of n and £ which enforce this mapping are 

fi = 2p-»/» and t=^r 3 . (2) 

When p* increases along the GCM freezing line, the BCC 
phase becomes eventually stable. The same effect occurs 
with the inverse-power and Yukawa potentials when 1 /n 
and I change with the freezing density according to Eq.|2 
In this respect, 1/n and £ play a role analogous to that of 
an effective temperature. For instance, dilute solutions of 
charged colloids with counterions , which constitute a 
practical realization of vi(r), show this kind of behavior 
as a function of the Debye screening length. 



The phase diagram of the GCM shows some resem- 
blance also with the phase behavior of star-polymer so- 
lutions ||. At low densities, the arm number / plays 
the role of an effective inverse temperature (for two rea- 
sons: the strength of the potential increases with / and 
its range increases with f^ 1 ). For not too huge values 
of /, the Yukawa repulsion yields a phase-stability sce- 
nario that is very similar to that of the GCM. Only for 
star-polymer packing fractions larger than ~ 0.7 (which 
is where the nearest-neighbor distance in a BCC solid is 
~ a) will the peculiarities of the short-distance repulsion 
make a difference, stabilizing other cubic phases that are 
likely unstable in the GCM. 

In this paper, we discussed the phase diagram of 
the GCM that was redrawn using current best-quality 
numerical-simulation tools. We predicted the existence 
of a narrow range of temperatures within which the se- 
quence of stable phases exhibited by the GCM upon 
isothermal compression is fluid, BCC, FCC, BCC again, 
and finally fluid again. We also rationalized these find- 
ings in terms of the properties of other softly repulsive 
potentials, with an emphasis on the phase behavior of 
star-polymer solutions. 

We acknowledge some useful discussions with Gian- 
piero Malescio. 
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